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Abstract 

The interaction between the Earth's magnetic field and the solar wind plasma re- 
sults in a natural plasma confinement system which stores energy. Dissipation of 
this energy through Joule heating in the ionosphere can be studied via the Auro- 
ral Electrojet (AE) index. The apparent broken power law form of the frequency 
spectrum of this index has motivated investigation of whether it can be described 
as fractal coloured noise. One frequently-applied test for self-affinity is to demon- 
strate linear scaling of the logarithm of the structure function of a time series with 
the logarithm of the dilation factor A. We point out that, while this is conclusive 
when applied to signals that are self-affine over many decades in A, such as Brow- 
nian motion, the slope deviates from exact linearity and the conclusions become 
ambiguous when the test is used over shorter ranges of A. We demonstrate that non 
self-affine time series made up of random pulses can show near-linear scaling over 
a finite dynamic range such that they could be misinterpreted as being self-affine. 
In particular we show that pulses with functional forms such as those identified by 
Weimer within the AL index, from which AE is partly derived, will exhibit nearly 
linear scaling over ranges similar to those previously shown for AE and AL. The 
value of the slope, related to the Hurst exponent for a self-affine fractal, seems to 
be a more robust discriminator for fractality, if other information is available. 



1 INTRODUCTION 

The characterisation of global energy transport in the coupled solar wind- 
magnetosphere-ionosphere system is a fundamental problem in space plasma 
physics EJ. Solar wind energy is transferred to, stored by, and ultimately re- 
leased from the magnetosphere by a range of mechanisms, in which substorms 
play a key role. Most investigations of the substorm problem have focused on 
single substorms or small groups of similar events, analogous to the study of 
individual earthquakes in seismology. 
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A complementary approach is to analyse inputs to and outputs from the 
system in an attempt to constrain the range of possible physics occurring in 
the magnetospheric "black box" (c.f. analogous approaches in climatology and 
seismology □). Reviews of the significant progress made so far in applying the 
methods of low dimensional chaos to the magnetosphere are given by Klimas 
et alu and SharmaQ; while more recent investigations into whether OLnot the 
"black box" can be treated as a self-organised critical (SOC) systemU are re- 
viewed by Watkins et alB, Chapman and WatkinsQ and Consolini and Changla. 
One mechanism for dissipation of magnetospheric energy is through Joule heat- 
ing in the ionosphere's auroral electrojets. This process can be studied via the 
auroral electrojet (AE) index, which is a means of estimating the electrojet cur- 
rent. The Joule energy dissipated depends upon both this and the ionospheric 
conductivity. AE is available at 1-minute resolution. Tsurutani et al.B showed 
this to have a "broken power law" frequency spectrum. The high frequencies 
approximately follow f~ 2 while the lower frequencies are / _1 with a break at 
about 1/5 h -1 . Power law frequency spectra are common in nature and can 
have several causes U such as Kolmogorov turbulence or the bifurcation route 
to chaos. They are thus in themselves not sufficient to completely constrain 
simple models. A parallel effort to studies of the power spectrum has been 
the search for Iqisl dimensionality, initially through the Grassberger-Procaccia 
(GP) algorithmoa. However, as noted by Osborne and ProvenzalelilJ, a low and 
fractional GP dimension is not uniquely a signaturxLof low dimensional chaos. 
It is also compatible with self-affine coloured noiset£l or SOCO. In view of the 
fact that AE is known a priori to be the output of a complex system, Takalo 
and Timonen, in an important series of papers 12-16 , investigated whether the 
dynamics of magnetospheric and auroral indices were better encapsulated by. 
stochastic "coloured noise" rather than by chaos. One test applied to A_EE3 
was for self affinity - a property of both coloured noise and chaos. A partic- 
ularly important technique for identifying self-affinity in the work of Takalo 
and Timonen 12-16 was the use of the second order structure function 5 2 (al- 
though other methods have also been applied to this problem 17-19 ). In this 
paper, by constructing a simple example, we illustrate that S2 alone cannot re- 
liably distinguish exponential autocorrelation from intrinsic self-affinity in the 
short timescale part of-the AE signal, which has been linked to the substorm 
"unloading" timescale By considering how 52 is related to other measures 
of self-affinity we address the question of what additional knowledge may be 
required to make S2 more useful. 
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Figure 1: Power spectrum of model Brownian motion (H = 0.5). 

2 SELF- AFFINITY (H = 0.5) IN BROWNIAN MOTION 

There are two kinds of fractal: self-similar and self-affme 0. They are dis- 
tinguished by whether the rescaling necessary to produce the original object 
is isotropic (self-similar) or anisotropic (sclf-affine). In the case of a random 
fractal such as a time series X(tL one is testing for statistical rather than exact 
self- affinity, so the test appliedE-3 uses the second order structure function Ed'Ej 
52(A), defined by 

Sa(A) =< (X(t + XAt) -X(t) f > (1) 
where < ... > denotes an average over time t. For a self afhne curve X(t), 

S 2 (X) ~ X 2H S 2 (1) (2) 

where H is the Hurst exponent IQ < H < 1 for a self-affine fractal) and 
52(1) =< (X(t + At) - X(t)) 2 >E3. This results in linear dependence (with 
slope H) of log [S'(A)/5'(1)] 1 / 2 on the logarithm of the_dilation factor A. We 
note that not only is it not necessary for A to be small li-3 but that self- affinity 
in fact implies that the above holds for all scales A. The time stationarity 
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Figure 2: Autocovariance (top panel) and scaling plot (bottom panel) for model Brownian 
motion (H = 0.5). 



assumption implicit in equation (1)0 allows us to use the definition of the 
normalised autocorrelation function ACF(XAt): 

ACF(XAt) . < ffl^gaa > (3) 



to rewrite S , 2 (A)/S , 2(1) in terms of the ACF 

S 2 (A) _ (1 - ACF(AAt)) 
l^OT ~ (1 - ACF (At)) ' 



(4) 



Alternatively one may form the numerator and denominator of (3) from the 
time-averaged, time-laggecL_products of the series SX = X(t) — Xlsee equation 
(1) of Takalo and TimonenEfl). We follow engineering conventiontiJ in referring 
to equation (3) with SX replacing X as the normalised autocovariance (ACV) . 
Equation (4) holds with ACVlXAt) replacing ACF(XAt), so either can be 
used as a tesj— for self-affinity 0-3. In numerical work we will follow Takalo 
and Timonenll3 in using the ACV. It is calculated for a discrete series (Xi\ 
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i = 1, N with mean X) by 

ACVW ^^-^g-V . (5) 

A classic example of a process which is both self-affine and fractal is Brow- 
nian motion. Figure 1 shows a representative power spectral estimate (unwin- 
dowed periodogram) for a time series of 131072 points of simple Brownian 
motion (H — 0.5). The well-known f~ 2 form is easily seen, limited only by 
the available dynamic range of the data. The upper panel of figure 2 shows 
the normalised autocovariance of the same time series. 

The lower panel of figure 2 shows log[S , 2(A)/S , 2(l)] 1 / 2 versus log A, where 
we calculate S2 using the normalised autocovariance from equation (5). The 
range in the plot of.-i|2 was chosen for ease of comparison with figure 4 of 
Takalo and TimonenOJ and our figure 6. The curves in both panels of figure 2 
are nearly straight lines. The value H — 1/2 can be read off from the slope of 
the line in the lower panel of figure 2. As expected, the structure function is 
an effective detector of its original intended target, a wide spectrum self-affine 
fractal signal. 



3 APPARENT SELF-AFFINE FRACTALITY (H = 0.5) IN EX- 
PONENTIALLY CORRELATED RANDOM PULSES 

The identification problem of self-affinity over a finite range begins to be ap- 
parent when one applies the structure function method to a series of random 
pulses. We first consider the case of random time series which have exponential 
autocorrelation function. Many physically-interesting random processes can be 
well approximated by an exponential AGFE3 As an exactly soluble example we 
note the simple "random telegraph" Ere 2 ]. This is a two level Poisson-switched 
process which switches between level £1 and level —F with constant probability 
1/r per unit time. This process hasLrEj an autocorrelation function: 

ACF(XAt) ~ e - 2 l AAt l/^ (6) 

which, by the Wiener-Khinchinc theorem, indicates a power spectrum of the 
form /~ 2 for high frequencies (/ ^> 1/r), but flat (/°) for low frequencies 
(/ < 1/r) & Bec ause e" 2 l AAt l/ T = 1 - 2|AAt|/r + 0(A 2 A 2 i 2 ) the scaling of 
log \/S2(A)/52(l) versus A will not only be linear (i.e. apparently self-affine) 
for XAt small compared with r/2 but will also give rise to a Hurst exponent 
value of 1/2 ii H is derived from the slope of the line (i.e. apparent fractality). 

Without knowing a priori that it is a 2-level, Poisson-switched system, ap- 
plication of S2 to a time series that was exponentially autocorrelated over time 
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could cause one to infer (erroneously) that the short lag behaviour correspond- 
ing to times AAi < r/2 was both self-affine and fractal. This serves to underline 
the point that self-affinity is an intrinsically wide bandwidth property, and that 
application of a wide-band test over the restricted range (XAt < r/2) makes 
it hard to distinguish certain types of randomness from self-affine fractality. 

4 APPARENT SELF- AFFINITY IN WEIMER PULSES. 

The relevance of the above observations to the AE time series becomes clearer 
when we consider that AE contains recurring "pulses" associated with magne- 
tosphcric substorms. Both the pulse shape and its recurrent properties could 
give rise to the observed scaling in AE. We first consider apparent scaling due 
to the pulse itself, and then examine the behaviour of a random series of such 
pulses. 

4-1 Restricted range self- affinity from a single Weimer pulse 

The pulse shape was studied by Weimer El in the AL index, one of the two 
indices from which AE is derived (AE — AU — AL). A random sample of 
55 substorms was divided into three classes based upon the peak AL value 
attained. For each class, the AL time series were superposed with respect to 
the substorm epoch, from which the average time series was then calculated. 
The three resultant average substorm profiles were shown to be well fitted 
by the functional form apte~ pt with both a and p increasing with increasing 
peak AL. This functional form is the solution of an ordinary second-order 
differential equation that was argued to describe the evolution of the electric 
field and currents in the substorm current wedge. The ionospheric part of the 
substorm current wedge is a westward current that the AL index was designed 
to measure. 

We now show that this shape causes apparent scaling in S2 at small values 
of XAt in the case of a single, isolated Wcimcr pulse. We take a = 1 without 
loss of generality. The numerator (ACF*) of equation (3) becomes: 

AC F* (XAt) =< pte~ pt p(t + AAi)e- p(t+AAt) > . (7) 

By starting with the identity 

e~ 2pt dt, (8) 

we may evaluate averages such as (7) by differentiation with respect to p. We 
find 

ACF*(XAt) = i(l + P XAt)e- pXAt (9) 
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and so using the denominator of (3) to normalise the ACF we have 



ACF(XAt) = (1 + P XAt) e - pXAt 



(10) 



Expanding the normalised ACF as a Taylor series gives 



ACF(XAt) = 1 - -p 2 X 2 At 2 



+ 0(A 3 ) 



(11) 



which yields, on insertion into the right hand side of equation (4), a scaling 
of 5 I 2(A)/S , 2(1) with A 2 , for XAt small compared with 1/p (observed to be 
« 30 minutes). This implies linear behaviour when the logarithm of either 
S I 2(A)/5 I 2(1) or its square root is plotted against log A. Hence the pulse then 
appears self affine over this range with H = 1. 

4-2 Restricted range self- affinity from random Weimer pulse train 

Now let us investigate the scaling properties of a sequence of such pulses, as 
might occur in the AE time series when measured, for .example, over the 100 
days (144000 points) studied by Takalo and TimonenEJ. As in the random 
telegraph we chose a random sequence of pulses, specialised here to a repre- 
sentative example of the Weimer pulse shape. Each pulse was of form apte~ pt 
where p — 1/30 minutes" 1 , a = 1, and the sampling interval was 1 minute 
for 131136 points. The inter-pulse intervals were drawn from an exponential 
distribution with e- folding time 300 minutes 24 . The above model is not meant 
to provide an exhaustive model for the AE time series, but the pulse is a 
knownE3 component of the AL (and thus AE) signal and so its contribution 
to the apparent self-affinity of AE must be investigated. 

Figure 3 shows a spectrum estimate for the model time series. The spec- 
trum has the characteristics of the exponentially autocorrclatcd random tele- 
graph with a breakpoint at around 1/p between f° for / < 1/p and / _4 for 
/ 3> 1/p. The time series gives rise to an autocovariance function with a steep 
(quadratic) slope at small lags XAt < 30 min(= 1/p) (see figure 4) character- 
istic of the pulse shape. The associated structure function has slope « 1 for 
XAt less than 10 minutes, and progressively less than 1 as XAt increases, such 
that it appears nearly linear over two decades in XAt (figure 5). 

Again this near-linearity, used alone without other information on a nat- 
ural signal of necessarily restricted dynamic range, could lead one to infer 
self-affine properties (or indeed chaotic ones) in a signal that is not self-affine. 
The addition of randomness to the single-pulse behaviour described in section 
4.1 has given rise to a Hurst exponent less than 1, when measured over the 
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Figure 3: Representative example of a spectrum from a random Weimer pulse train. 



whole of the range 1 < XAt < 100. We believe there to be competition be- 
tween the effects of randomness (e.g. H = 0.5 in the random telegraph) and 
the integer value of H = 1 associated with individual pulses. 

5 AE RE-EXAMINED, DISCUSSION AND CONCLUSIONS 

We now consider the scaling properties of the measured AE time series in the 
light of the previous examples. The top panel of figure 6 shows the autocovari- 
ance of the first 100 days of AE for 1983_and may be compared with the top 
panel of figure 4 of Takalo and Timoneno. Again, the steep (exponential) fall 
of the ACV results in a near linear slope for small AAi, and a slow decrease 
in the slope as larger and larger ranges of AAi is considered . Importantly, 
however, the slope is, always less than 1 (J. Takalo, Private communication, 
1999). Overall it resembles near-linear scaling in the structure function for 
most of the first two decades of A (bottom panel of figure 6), and (plotted in 
the the middle_Danel of figure 4 of Takalo and Timonentj), was cited by Takalo 
and TimonenE-3 as a key piece of evidence for self-affinity in AE. They also 
noted the resemblance of the AE autocovariance function to an exponential 
and proposed that the autocorrelation time of AE be defined as the lag for 
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Figure 4: Autocovariance of random Weimer pulse train 



100F 




1 10 100 1000 

A (minutes) 

Figure 5: Scaling plot for random Weimer pulse train. 

timescales longer than that over which the autocovariance ceased to be expo- 
nential. Inspection of figures 4, 5 and 6 lead us to conclude, however, that, 
unlike the ideal case of Brownian motion, neither the curve of S2 for AE, nor 
that of the simplified random Weimer pulse train are straight over the range 
XAt = 1 to A = 120. We remark that, insofar as the ACF of AE is exponential 
for small XAt, there must eventually be a departure from near-linearity in the 
structure function as XAt increases, unless the range over which the exponen- 
tial behaviour is seen is so small that a straight line would be just as good an 
approximation as the exponential. 

In addition, both AE and the model Weimer pulse train of section 4.2 
give a fractional H value when taken over the whole range from 10 to 100. 
Without a priori additional knowledge, we might equally well have concluded 
that the random pulse train was self-afline over the range XAt < 100, but by 
construction we know this is not so. 

Our model was deliberately simplified. In the natural AE time series, the 
extended tail of the ACV is expected to reflect the solar wind-driven component 
(also present in AU and AL), -which our simulation neglected. As originally 
conjectured by Tsurutani et ala the solar wind driver is probably the origin of 
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Figure 6: Autocovariance and scaling plot for 100 days of AE, starting on 1st January 1983 

the "1//" part of the AU/AL/AE spectrum!!!]. 

We may summarise our findings as follows. By construction of an explicit 
counter-example we have shown that near-linear scaling of S2 over about two 
decades is not in itself sufficient to show self-affinity. We have also given 
analytic and numerical evidence that non-fractal random series can produce 
non-integer Hurst exponents over limited dynamic ranges. We thus infer that 
self-affinity in the range to 100 minutes for AE has not been and could not 
be proved by the use of S2 alone. 

One may reasonably point out that several other methods have been used 
to provide evidence of self-affinity in geomagnetic indices; both in the papers 
of Takalo et al. and those of other workers 17-19 . One may thus enquire as to 
what kind of additional knowledge or analysis techniques would be necessary 
for considering the results of the structure function method fruitful? Based on 
what we have found, we suggest that answer is at least threefold. 

1) Be aware that many tests for fractality are actually designed assuming 
a fractal signal: A test based on the assumption of fractality can disprove 
fractality but cannot prove it. The methods for measuring fractal dimension 
that we are aware of assume self- affinity in their design i.e. they typically 
examine the scaling behaviour of a signal. Only if they find no evidence of 
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scaling at all is there no ambiguity. 

2) Use more than one test: Several tests are better than one because 
different methods are sensitive to different non-fractal effects. Thus use of 
several tests means that a series with non-fractal aspects is less likely to be 
misinterpreted. Most of the methods for measuring fractal dimension which 
have been applied to geomagnetic data are of one of two basic types. The 
first type of method basically estimates the dimension of a fractal curve by 
examining how the average value of short lengths of curve 

Si =< X(t + XAt) -X(t) > (12) 

scales with the ruler length A (in units of the sampling interval At). Such 
methods have been applied by VorosllH to magnetometer data, and mare re- 
cently to geomagnetic and solar wind quantities by Price and Newmanl23, who 
used the related, cumulative "R/S" analysis. 

The second type studies the positive definite second order function 

S-2 =< (X(t + XAt) — X(t)) 2 > (13) 

and returns the same information!^ as the ACF when estimated on a stationary 
signal (see section 2). For this reason it is thus also formally related to the 
power spectrum via the Wiener-Khinchine theorem. S2, the ACF and the 
power spectrum have all been extensively investigated for the AE, DSt and 
related indices by Takalo et al 12 ~ 16 . The meaning of this family of techniques 
can be understood as studying the behaviour of the histogram of variance 
of the signal (or the power spectrum) with increasing time dilation XAt (or 
frequency); depending on whether one is dilating in time (in the case of S% 
and the ACF) or frequency (in the case of the power spectrum). We caution 
that time lag in the ACF or in S2 is not trivially l/(the Fourier frequency) 
because any frequency in a Fourier transform has contributions from multiple 
ACF lags and vice versa (see Bendat and PiersolEil, pages 120-122). In the 
case of a simple fractal, the dimension (and Hurst exponent H) estimated 
from such methods should theoretically be the same as frorrLjS'i, although in 
practice the errors of the two methods need not be the samelia. If they differ 
substantially, this may be a pointer that the time series is not intrinsically a 
wideband fractal, and that one of S\ or S2 is more sensitive to this. 

An example of how additional tests for fractality have supplied new knowl- 
edge is in the continuing study of the AE index. This has been known since 
the work of Tsurutani et ala to have a "1//" low frequency and "l// 2 " high 
frequency power spectral density. Acting only on information from the power 
spectrum or other SVtype methods, one might thus infer that AE is a bi-affine 
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quantity 12-16 , i.e. it has two_separate scaling regions and a break hetpzeen 
them. In contrast, ConsoliniEj has studied the "burst distributions" ESEj of 
AE. These are the histograms of intervals between threshold crossing times 
(burst and inter-burst lifetimes) and of areas above threshold between c:pss= 
ings (burst sizes). The lifetime distributions are an Si -type measurement Eo 
and were found to have (exponentially rolled-off) scaling with a single slope 
over a very wide range, interrupted only by a non-scaling component at about 
2 hours. The apparently paradoxical observation of bi-affine behaviour in S2 
and "contaminated" mono-affine behaviour in Si has been addressed in two 
different ways. One has been to introduce models which—have the required 
properties in both Si and S2, such as forest fire models c§ or coupled map 
latticesc3 driven by wideband solar-wind like signals. The other, informed by 
the fact (section 4 and 5 above) that the high frequency f~ 2 part of a power 
spectrum jaeed not arise from a fractal aspect of the time series, has been to 
postulate that the AE series is in fact a hybrid time series with a fractal 
element arising from the solar-wind driven ionospheric current systems and a 
non-fractal part arising from energy storage and release in the magnetosphere 
(substorms). This was supported E3 by the observation that the scaling in AE 
(and AU/AL) burst-lifetimes is the same as that seen in the solar wind (see 
also Freeman et alEll) while the non-scaling component was seen only in the 
magnetospheric quantities such as AU and j4Xcj. 

There have been exceptions to the use of S? or Si type techniques in the 
geomagnetic context. We are grateful to an anonymous referee for reminding 
uSpaf the results of a multifractal analysis of the AE index by Consolini et 
al.H. These results must imply some constraints on possible models describing 
the variability of auroral currents. However, in the same way that AELwhen 
measured over 1 year by a method of type Si is essentially fractal cJ, and 
required the use of several years' measurements for the "bump" -like feature in 
the otherwise scaling Si to become apparent E3, it seems tp_.us that one might 
expect a multifractal analysis of less than 2 months of AE'Ej to give a good fit 
to a p-model of turbulence because the solar wind driver is also well fitted by 
this particular turbulence model E3. We believe that a study on a much longer 
series of AE would be required to exclude even our own toy model of random 
differentiable Weimer pulse trains, when superposed on the multifractal solar 
wind background. We note that, independently, a recent multifractal study of 
geomagnetic data from Thule, Alaska has excluded the biaffine coloured noise 
model E-3 for that dataset. 

3) Remember that Nature does not have to be purely fractal any more than it 
has to be non-fractal: Many types of natural signal have both fractal and non- 
fractal components. In consequence, when using methods to examine fractality, 
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one should be aware that it is possible to find something between the extremes 
of wideband fractality and none at all, as discussed in point (2) above for the 
case of AE. Another example is to imagine looking out of one's tea-room 
window at a tree through a regularly spaced window blind. The distinguishing 
of the fractal tree and the periodic blind is a task that the human eye and 
brain perform routinely, and which a Fourier transform can also do because 
it can resolve the blind spacing as a spatial frequency. A "random blind" 
appearing at Poisson-switched intervals would be much more of a problem for 
an FFT, and would be analogous to the pulses of section 3 and 4. The user 
thus needs to determine how much the presence of a "contaminating" signal 
or signals in the fractal time series may affect their interpretation, at which 
point the question may become as much physical as mathematical. This is 
currently an admittedly very difficult task because of the sparsity of literature 
on such hybrid time series, and is one which we plan to examine in more detail 
in future papers. 
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